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Abstract. A 34'x34' field centred on the spiral galaxy 
M 81 has been searched for emission-line objects using 
the prime focus wide field camera (WFC) of the 2.54 m 
Isaac Newton Telescope (La Palma, Spain). A total of 171 
candidate planetary nebulae (PNe) are found, 54 of which 
are in common with the ones detected by Jacoby et al. 
(1989). The behaviour of PNe excitation as a function 
of galactocentric distance is examined, and no significant 
variations are found. The PNe luminosity function is built 
for the disk and bulge of M 81, separately. A distance 
modulus of 27.92 ± 0.23 mag is found for disk PNe, in 
good agreement with previous distance measurements for 
M 81 (Jacoby et al. 1989; Huterer et al. 1995). 
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1. Introduction 

Basic information about planetary nebulae (PNe) in ex- 
ternal galaxies has been summarized by Jacoby (1997). 
Many hundreds of PNe were discovered in spiral galax- 
ies (mainly in M 31), in irregular (LMC and SMC) and 
elliptical galaxies up to the distance of the Virgo cluster. 
More recently, 131 candidate PNe have been discovered by 
Magrini et al. (2000, 2001) in M 33. PNe are important 
objects in the study of stellar population of intermediate 
age in galaxies of different morphological types and in dif- 
ferent chemical environments, and in the assessment of the 
kinematical properties of all morphological components of 
galaxies (disks, bulges, haloes). 

M 81 (NGC 3031), a nearby galaxy outside the Local 
Group, belongs to the group of the most massive spiral 
galaxies (SAab) being also the nearest LINER galaxy. To- 
gether with M 82 (NGC 3034), NGC 3077 and various 
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dwarf galaxies, it forms a quite interesting multiple sys- 
tem. The central area of 4'x4' of M 81 was searched for 
PNe by Jacoby et al. (1989), who discovered 185 candi- 
date PNe. In the present paper we address the search of 
PNe in a much more extended region, a factor of 70 times 
larger than the already surveyed area, albeit with a 0.7 
mag lower detection limit in the central crowded regions. 
Our survey covers the whole optical extent of the galaxy 
and its surroundings. In Sect. 2 observations and data re- 
duction are presented. Sect. 3 contains the data analysis. 
In Sect. 4 we introduce the newly discovered candidate 
PNe, while in Sect. 5 we discuss the behaviour of their 
excitation in comparison with other galaxies. In Sect. 6 
we build the PNe luminosity function in order to evaluate 
the completeness of our sample and to use it as a standard 
candle for extragalactic distance determination, thus de- 
riving the corresponding distance modulus of M 81. The 
summary and conclusions are presented in Sect. 7. 

2. Observations and data reduction 

The M 81 galaxy was observed on December 14, 2000 and 
January 2, 2001 using the prime focus wide field cam- 
era (WFC) of the 2.54 m Isaac Newton Telescope (La 
Palma, Spain). An area of 34'x34', covering the whole 
galaxy (see Fig. 1), was surveyed with a panoramic de- 
tector consisting of four EEV CCD of 2048 x 4096 pixels 
each. The pixel size projects to 0".33 in the sky. Observa- 
tions were taken through three filters with the following 
central wavelengths and widths (FWHM): an [O III] fil- 
ter (5008/100 A), an Ha+[N n] filter (6568/95 A) and a 
Stromgren Y filter (5550/300 A), the latter used as an off- 
band filter for identifying the emission- line sources. The 
total exposure through the [O in] filter was of 6600 s split 
into two sub-exposures of 40 minutes and one of 30 min- 
utes. The Ha+[N n] image had a total exposure time of 
1800 s (three 600 s frames) whereas the Stromgren Y im- 
age was exposed 1600 s (four 400 s exposures). The centre 
of the telescope pointing was: 09h 55m 33.2s, +69d 03m 
55s (J2000.0). The seeing was approximately 1".5 in the 
[O in] and Strom. Y frames, and 1" in the Ha+[N n] 
frames. The CCD frames were de-biased and flat-fielded 
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using IRAF. All images were aligned to a reference image, 
correcting for geometrical distortion. Images in the same 
filter were then averaged out to remove cosmic rays and 
to improve the signal to noise ratio. Sky background was 
subtracted using the external regions of our frames, where 
the optical emission of the galaxy is negligible. 

3. Data analysis 

Emission-line objects were identified by subtracting the 
stellar continuum of M 81 using the off-band, Strom. Y 
frames from the [O in] and Ha+[N n] frames, as described 
in Magrini et al. (2000). This procedure removes the stel- 
lar background allowing for identification of emission-line 
sources. 

Following Jacoby et al. (1989), who tested the dif- 
ferences between crowded field photometry versus aper- 
ture photometry for the M 81 bulge PNe, we adopted 
the aperture photometry technique. [O in] and Ha+[N n] 
line fluxes were then measured using APPHOT in the 
continuum-subtracted images. Errors were estimated con- 
sidering both photometric errors, given by Poissonian 
statistics on the background and on the sources, and the 
detector noise. Photometric errors vary between few per 
cent for the brightest objects in the [O in] and Ha+ [N ii] 
images up to about 20% for the fainter ones in the [O in] 
frames, whereas they can reach 50-60% for the faintest 
candidate PNe in the Ha+[N n] images, where the sig- 
nal to noise ratio is lower due to the shorter exposures. 
Instrumental fluxes were transformed into physical fluxes 
using the standard PN calibrator PNG 205.8-26.7, whose 
accurate fluxes are given in Dopita & Hua (1997). 

Emission-line fluxes were then corrected for the in- 
terstellar extinction. Following Jacoby et al. (1989), we 
consider only the foreground Galactic extinction towards 
M 81, which is quoted to be Ay — 0.3 mag, and corre- 
sponds to E(B — V) ~ 0.1 mag (Kaufman et al. 1987). 
We have adopted this value to correct the [O m] and 
Ha+[N n] fluxes according to Seaton (1979). 

Coordinates of the emission-line sources in our images 
were derived with a multi-step procedure. Firstly, approx- 
imately 30 stars for each [O in] CCD field were identified 
on the Digitized Sky Survey []. A first astrometric solution 
was computed with these stars, and with this solution co- 
ordinates for about 400 stars found with DAOFIND in 
the field were subsequently obtained. Their coordinates 
were then replaced with the nearest APM-POSS1 coordi- 
nates and a new astrometric solution was computed. The 
procedure was iterated to reach a final accuracy of ap- 
proximately 0."5 r.m.s. 



1 Based on photographic data of the National Geographic 
Society Palomar Observatory Sky Survey (NGS-POSS) ob- 
tained using the Oschin Telescope on Palomar Mountain. 



4. Candidate planetary nebulae 

To identify PNe, we have adopted the following criteria: 

i) the object should appear both in the [O III] and 
Ha+[N n] frames but not in the continuum frame. 
A negligible continuum is in fact expected from a PN 
at these wavelengths; 

ii) it should not be spatially resolved, i.e. it should have a 
stellar point spread function. At the distance to M 81 
(3.50±0.40 Mpc, Jacoby et al. 1989), 1" corresponds 
to approximately 17 pc whereas the largest galactic 
PNe have sizes of ~4 pc (cf. Peimbert 1990; Corradi 
et al. 1997). Therefore, all PNe belonging to M 81 are 
expected to be unresolved in our images, and any ex- 
tended emission region is considered to be an H II re- 
gion or a SN remnant. 

One hundred and seventy-one objects fulfilling criteria i) 
and ii) have been found. They are listed in Table 1 to- 
gether with their observed [O in] and Ha+[N n] fluxes 
and coordinates; they are indicated in Fig. 1 by "+" signs. 
Fifty four candidates coincide with the objects discovered 
by Jacoby at al. (1989) in the bulge of M 81; their coordi- 
nates and fluxes agree with those in Jacoby at al. (1989) 
within the errors. In column 6 of Tab. 1 the identification 
number in Jacoby et al. (1989) is reported. 

Magrini et al. (2000) have shown that the 
R=I([0 Ill])/J(Ha+[N ii]) flux ratio can be used to 
distinguish PNe from other emission-line objects in a 
statistical sense, but not for individual objects. In fact, 
the analysis of the catalogue of Galactic PNe by Acker 
et al. (1992) shows that R spans a wide range, roughly 
from ([O in] undetected) to 10, and that ~75% of 
PNe have R > 1. H II regions are instead generally of 
lower excitation (R < 1), but this is not always true, 
so that misclassification between PNe and compact H II 
regions (which appear as point-like at the distance of 
M 81) can occur. Only detailed spectroscopy would allow 
a final distinction between the two classes of objects. 
Nevertheless, using the R flux ratio it is at least possible 
to estimate the contamination of H II regions in our 
sample in a statistical way, as done by Magrini et al. 
(2000). For the 135 objects in Table 1 in which both 
the [O in] and Ha+[N n] fluxes have been measured, 
we find 28% with R < 1, as compared to 25% for the 
Galactic PNe (808 PNe whose fluxes are quoted in Acker 
et al. 1992). Assuming that M 81 and the Milky Way 
contain the same percentage of low excitation PNe, then 
we expect only a few H II regions (approximatively 5) 
contaminating our sample of candidate PNe. 

5. Radial distribution of excitation and 
comparison with other galaxies 

The excitation class of a PN is defined in terms of flux 
ratios, primarily of [O in] to Ha or H/3 but also consid- 
ering other nebular lines (Feast 1968). These ratios are 
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Fig. 1. The Palomar Atlas image of M 81. The field of this frame is approximately 30'x30'. The field of view of the 
INT+WFC is slightly larger (34'x34'), but we show here only the area in which we have found PNe. North is at the 
top, East to the left. The "+" symbols indicate the location of the candidate PNe. 



good indicators of the temperature of the central star 
but also depend on the properties of the nebula like its 
geometry, electron density and temperature (and there- 
fore the abundance of the important coolant O/H) and 
on the optical depth in the Lyman continuum. The flux 
ratio R = h ^+[n\i\ thus provides us with an indication 
of the excitation of the candidate PNe in M 81, as already 
anticipated in the previous section. 



In the Galaxy, Vorontsov-Vel'yaminov et al. (1975) 
presented an absolute spectrophotometry of 47 PNe in the 
direction of the Galactic centre. They found that their R 
ratios differ significantly from those of disk PNe, conclud- 
ing that there are physical differences between Galactic 
bulge and disk PNe. Webster (1975, 1976) compared the 
distribution of excitation of the Galactic bulge PNe with 
those in the Magellanic Clouds and proposed the presence 
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Fig. 2. The flux ratio R = u a +[Nii ] as a function of 
the galactocentric distance for the Milky Way (data from 
Acker et al. 1992, upper box), M 81 (this work, middle 
box), and M 33 (Magrini et al. 2000, Magrini et al. 2001, 
lower box). 



of more massive progenitors in the case of LMC, where 
there are more high excitation PNe than low excitation 
ones, compared to our Galaxy and the SMC. In spiral 
galaxies, information on the excitation of PNe thorough 
the bulge and the disk is available only for M 81 (this pa- 
per), M 33 (Magrini et al. 2000) and the Galaxy (Acker 
et al. 1992). Using the information contained in the Stras- 
bourg ESO Catalogues, the galactocentric distances of 427 
Galactic PNe with quoted Ha, [N n] and [O in] fluxes 
were computed; the distances from the Sun of the PNe 
come from Cahn et al. (1992), and the adopted galacto- 
centric distance of the Sun is 8.5±1.1 kpc (Allen, 2000). 

Fig. H shows R as a function of the distance from 
the centre of the galaxy in the three cases. PNe in M 33 
and the Milky Way are approximately homogeneously dis- 
tributed in terms of excitation. The apparent progressive 
reduction of high excitation PNe with galactocentric dis- 
tance in M 81 (Fig. 2) is removed if the mean R, computed 
in bins of 2 kpc, is plotted (Fig. 3). 

Therefore, and contrary to Vorontsov-Vel'yaminov et 
al. (1975), no significant difference between the excitation 
of bulge and disk PNe is found in the Galaxy nor in M 81 
(Fig. I). 



Fig. 3. As in Fig. 2, but plotting the mean value of R, 
averaged on bins of 2 kpc and up to the distance of 12 
kpc. The error bars indicate the standard deviation from 
the mean. 

6. Planetary Nebulae Luminosity Function 

[O III] fluxes for the candidate PNe in M 81 were con- 
verted into equivalent V-band magnitudes following Ja- 
coby (1989): 

m [0 in] = -2.51ogF [0 up] - 13.74. (1) 
The Planetary Nebulae luminosity function (PNLF) 
was then built for three different samples: all 171 PNe 
of M 81, those in the disk (86 objects) and those in the 
bulge (85). They are presented in Fig. |[ The Eddington 
formula (Eddington 1913) was applied to correct for the 
effects of observational errors. The resulting luminosity 
function within the completeness limit was then fit to the 
"universal" PNLF 

%)ae aw(m »(l-e 3(ra '- m) ). (2) 
In Eq. (2) m is the observed [O in] magnitude from Eq. 
(1), and m* is the apparent magnitude of the PNLF cutoff 
of M31 (Jacoby 1989). As we have discussed in Magrini et 
al. (2000), it is possible to estimate the limit of complete- 
ness of a sample of PNe by computing the S/N ratio for 
each PNe (putting the limit of completeness at S/N~10, 
Ciardullo et al. (1987)) or equivalently by evaluating the 
point at which observational PNLFs inflect downwards in- 
dicating incompleteness. Fig. || shows that the complete- 
ness limit changes for the different samples, indicating the 
strong dependence on the galaxian background emission, 
which increases towards the centre of the galaxy. Com- 
paring our data with the sample of Jacoby et al. (1989), 
we note in fact that we are losing approximately 70% of 
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Fig. 4. Histograms of R for M 81 (this work, upper box) 
and the Milky Way (lower box), for disk and bulge PNc 
separately. 



the PNe found by Jacoby and collaborators in the bulge 
of M 81. This is due to two effects: the use of a filter 
that is three times broader than the one used by Jacoby 
et al. (1989), and a smaller primary mirror (2.5 m vs. 
4 m) not compensated by longer exposure times (1.8 hr 
vs. 1.5 hr). These effects are particulary important in the 
central regions because of the high surface brightness of 
the parent galaxy. In these regions our completeness limit 
of 24.0 mag differs from that of the sample of Jacoby et 
al. (24.7 mag) and from the one we have for our data in 
the disk (24.7 mag). For this reason the distance modu- 
lus of M 81 is computed using only disk PNe, where we 
consider to have a complete sample. The error on the dis- 
tance modulus is obtained by combining the error asso- 
ciated with our best least mean square fits (±0.04 mag) 
in quadrature with those associated with the photomet- 
ric zero point (0.03 mag) and with the adopted extinction 
(0.2 mag from Kaufman et al. 1987). We have also to con- 
sider two systematic errors which come from the uncer- 
tain definition of the empirical PNLF (0.05 mag) and the 
distance of the calibration galaxy M31 (0.10 mag). These 
errors lead to a total error of ±0.23 mag. The resulting dis- 
tance modulus to M 81 is 27.92 ±0.23 mag, corresponding 
to 3.84 ± 0.41 Mpc. This distance is in excellent agree- 
ment with the value derived from the Cepheids method 
of 27.97 ± 0.14 mag (Huterer at al. 1995) and also in fair 
agreement with the value derived by Jacoby et al. (1989) 
using the luminosity function for their bulge PNe candi- 
dates (27.72 ±0.25 mag). 



Fig. 5. Mean square fits for the luminosity functions of 
PNe in M 81: i) Complete sample (171 candidate PNe); 
ii) Disk (86 candidate PNe) ; iii) Bulge and central regions 
(85 candidate PNe). Black points indicate candidate PNe 
within the completeness limits. 

7. Summary and conclusions 

We imaged a 34'x34' area around the spiral galaxy M 81 
through [O in] and H«+ [N n] narrowband filters. In total, 
171 candidate PNe are found, 117 of which are new while 
54 coincide with those discovered by Jacoby et al. (1989) 
in the central 4'x4' bulge. 

The behaviour of PNe excitation across the galaxy was 
examined, finding no evidence for substantial differences 
in excitation between bulge and disk PNe, nor for vari- 
ations along the galaxian disk. The same result applies 
to the Milky Way and M 33, and is contrary to previous 
suggestions from Vorontsov-Vel'yaminov et al. (1975). 

The PNLF of the candidate PNe in the disk of M 81 
provides a distance determination of 3.84 ± 0.41 Mpc. 

Acknowledgements. We are grateful to Peter Sorensen for tak- 
ing the observations of M 81 in service mode. 
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Table 1. PN candidates found in M81. Fluxes in Ha+[N n] and in [O in] are in units of 10~ 16 erg cm~ 2 s _1 and are 
corrected for reddening (see text). 
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